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We study charge ordered solutions for fermions on a square lattice interacting with dynamic 
antiferromagnetic fluctuations. Our approach is based on real space Eliashberg equations which 
are solved self-consistently. We first show that the antiferromagnetic fluctuations can induce arc 
features in the spectral functions, as spectral weight is suppressed at the hot spots; however, no 
real pseudogap is generated. At low temperature spontaneous charge order with a d-form factor can 
be stabilized for certain parameters. As long as the interacting Fermi surfaces possesses hot spots, 
the ordering wave vector corresponds to the diagonal connection of the hot spots, similar to the 
non-self-consistent case. Tendencies towards observed axial order only appear in situations without 
hot spots. 

PACS numbers: 71.10.Fd,74.72.-h,71.27.+a 


I. INTRODUCTION 

Signatures of charge order in the copper-oxide based 
superconductors have attracted a lot of recent attrac¬ 
tion. By now charge ordered states have become an es¬ 
sential ingredient of the phase diagram for the different 
families.^ Whilst in La-based compounds charge and 
spin stripe order have a considerable historyP^^ it was 
only firmly established in recent years in other cuprate 
families. Early reports include scanning tunneling mi¬ 
croscopy (STM) studies,followed by bulk pro perty 
measurements by resonant X-ray scattering (REXS).^^^^^ 
The order has a finite correlation length in zero magnetic 
field and b ecomes long range ordered at high magnetic 
fields.11^112 This helps to understand previously puzzling 
quantum oscillation data.f^^H^ 

By now the charge order has been carefully character¬ 
ized experimentally. It is unconventional in the sense that 
it possess an incommensurate wave length A = 3 — 4a, 
where a is the Cu-Cu distance in the copper oxide planes 
and an internal form factor. In both STM and REXS 
the direction of the wave vector has been established 
as pointing along the axes in the Cu-Cu square lat¬ 
tice, Q = Qo(l, Its magnitude follows the 

Fermi surface, or more precisely, Fermi arc evolution; 
concretely this means that it decreases with increasing 
hole doping.^ This suggests that the order is connected 
to Fermi surface properties and ca n pos sibly be under¬ 
stood as a Fermi surface instability)In contrast, in 
La-based compounds the ordering wave vector shows the 
opposi te tre nd, and Q increases with increasing hole 
dopingRecent experiments have revealed further de¬ 
tails of the charge order in Bi- and Y-based compounds. 
Both STM and REXS data are well understood based on 
a (i-form factor as an internal structure for the order.l^^lt^ 

Many theories have addressed the microscopic origin 
of the charge order including its onset temperature, its 
wave vector, f orm factor and dependence on doping and 
magnetic field.One line of reasoning starts 



FIG. 1: (Color online) Schematic picture of Fermi surface with 
antiferromagnetic wave vector connecting hot spots K and 
examples of diagonal, Q — (3o(l, 1)? cind axial, Q — (3o(l,0), 
charge ordering wave vectors. 

with fermions interacting with antiferromagne tic sp in 
fluctuations peaked at wave vector K (see Fig.[^.l^^®E^ 
The hot spot model with linearized dispersion has a 
strong nesting instability with a diagonal wave vector 
Q = Qo{l,l) connecting the hot spots (see Fig. 

In models starting with an intact Fermi surface this in¬ 
stability is at least as strong as the one with an axial 
wave vector Q = QoCl^O), and it is therefore difficult to 
explain the experiments based on such a theory. A num¬ 
ber of recent works have explored situations where the 
non-interacting Fermi surface has been modified due to 
magnetic order and pocket formatiorP^®^ or pseudogap 
features.l^It has been argued that it is important to take 
such features into account, such that the dominant insta¬ 
bility can be alter ed. A lso fluctuation effects,l^^lfl3 strong 
correlation effects,®^ and starting points of fractional 
Fermi liquidP^ have been considered as a possibility to 
favor axial order over the diagonal one. 

The main idea of this work is to analyze whether sin¬ 
gle particle self-energy modification, due to a coupling to 
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dynamic spin fluctuations, and without appeal to more 
exotic states, is sufficient to establish charge order with 
properties as observed in experiments. Our analysis will 
be carried out with a semi-phenomenological model of 
fermions coupled to dynamic spin fluctuations, also called 
spin-fermion model in the literature. Such a model is ap¬ 
pealing as it naturally capt ures t he d-wave superconduc¬ 
tivity in a range of dopingsWe show that in the rel¬ 
evant doping regime for charge order, dynamic magnetic 
fluctuations with a finite correlation length ^ ^ 2a are a 
more realistic descrip tion f or the cuprates than including 
static magnetic order.E^EH Xo analyze the model, we use 
a real space formulation which allows us to probe charge 
order directly in the ordered phase. Our formulation also 
allows for simple extensions to include impurities or dis¬ 
order, spin order, coupling to phonons, superconducting 
order, and a magnetic field. The paper is structured as 
follows: In Sec. II we describe the details of the model 
and approach. In Sec. Ill we analyze ordered solutions 
in the static limit in order to connect to previous work. 
In Sec. IV we analyze the full dynamic model before con¬ 
cluding in Sec. V. 


II. MODEL DEFINITION AND REAL SPACE 
ELIASHBERG APPROACH 


We use a model with the action of the form,^^^^ 

S = ~ ^ ^ V^i,cr(^‘^n)G^27\Q(^‘^n)'0ji',cr(^‘^n) (1) 

(j ij,n 

- y V xin - rj,T - T')S{ri,T) ■ S{rj,T'), 

hj 

where 

“ 1 “ ^ij ^ (^) 

and 

r) = y ('r)cr“^ ,<^2 (t). (3) 

Here, is a fermionic field for site i and spin label a. 
The non-interacting dispersion in momentum space is 


Eh = —2ti[cOs(/Ci) + COs(/C2)] — 4^2 COs(/Ci) COs(/C2) 
—2ts[cos{2ki) + cos ( 2 / c 2 )]. 


where the copper-copper lattice spacing is set to a = 1. 
For the spin-fluctuation spectrum we assume,^^^^ 
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(4) 


To preserve lattice periodicity we define uokAo) = 2(2 — 
cos(ga, — Kx^i) — cos{qy — Ky^i)) + T^, and we usually 
take for simplicity K = T is directly related to 

the correlation length T = For a^^ = 0 this form is 


common in the literature such as in the well-known work 
of Millis et al^ and fits experimental data from neutron 
scattering well. In real space we have 

X{ri - rj,iujm) = T y] x(g, (5) 

* Q 

where Ng is the number of lattice sites. 

a. Real space equations - The basic equation for the 
Eliashberg approach is 

YjijiiuJn) — ^ ^ X('^2 ^ (6) 

m 

where G is the full Green’s function. We assume that 
no spin order occurs and omit spin labels on G and E. 
In order to compute E self-consistently we also need to 
solve the Dyson equation, 

GijiiuJjn) — tij 

When working in real space Gij{iuJm) can be obtained 
by matrix inversion for each a computation which 

can be easily parallelized. The calculation of E in Eq. 
can also be parallelized well. In the Green’s function G in 
Eq. 0, we have included a general symmetry breaking 
field, which can be introduced into the Hamiltonian by 

Hqo = (^) 

In order to find inhomogeneous solutions we have to ini¬ 
tialize the calculations with such a field. We assume that 
X in Eq. Q is fixed by the input parameters and not fur¬ 
ther renormalized in a self-consistent manner. We work 
on a two-dimensional lattice with Ng = Ni x Ni sites 
and periodic boundary conditions. 

b. Homogeneous case - In the homogeneous situa¬ 
tion G and E only depend on Vi — Vj. With a Eourier 
transform we can then write. 


^k{i^n) = x(fc - 9, - iLOm)G{q, iuim) (9) 


and the Dyson equation can be solved explicitly, 

Gk{i^m) = ^ I T~. 7- 

In the spectral representation we have 


( 10 ) 


“ J J -UJ1-UJ2 

X [^f(<^2) + ^b(—^ 1 )], 

where p^{q,uj) is the spectral function for the spin- 
fluctuation spectrum, and nF(a;), nB{oj) denote the 
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fermionic and bosonic distribution functions, respec¬ 
tively. Analytic continuation, iuj ^ uj ^ irj, yields the 
imaginary part of the self-energy 




-Sng'^ y] Jduj2 Px{k - q,u} - uJ2)PG{q,i^2) 


x[nF(w2) + nB(w2 - w)]. 


( 11 ) 


and the real part can be computed from the Kramers- 
Kronig transformation. From this the electronic spectral 
function pcik^uj) can be computed. For a numerical cal¬ 
culation it is favorable to use the real space form for the 
self-energy equation and switch to momentum space by 
Fast Fourier transforms. 

c. Charge order observables - For the order param¬ 
eter we focus on the inhomogeneous part of 


^ij — 


cr Q k 


^iQR (-^2) 


the following parameters: the overall weight factor a^, 
the inverse correlation length F = the spin fluctua¬ 
tion scale cjgf and the coefficient . We use parame¬ 
ters similar to the recent work of Mishra and Norman.l^ 
In addition we have the temperature T, where the lowest 
value reached is T = 0.02ti 70K). Moreover we have a 

coupling constant g. A technical parameter is the size of 
the real space lattice. We did most of our calculations for 
Ni =32. The limit for this is set by memory constraints. 


III. RESULTS FOR THE STATIC MODEL 


We first consider the situation where the spin- 
fluctuation spectrum is a static, which corresponds to 
the limit cjgf ^ oo and ^ 0 in Eq. This helps 
us to connect to previous result^^ and test our formal¬ 
ism and procedure. The equation for the self-energy then 
simplifies to. 


where r = Vi — Vj and R = {vi rj)/2. Given a 
spin-diagonal Green’s function Gij{iuJn) we can com¬ 
pute ^ 0. To ex¬ 

tract the inhomogeneous part we map it to coor¬ 
dinates R, r, A^j A(R, r) and then compute A(r) = 
A(R, r) — r). Then we find A^^(R, r) = 

A(R, r) — A(r). Without any charge order symmetry 
breaking A^^ is zero. 

We assume that the inhomogeneous part Aq(/c) can 
be expanded as 

^Q{k) = Y.an{Q)Uk), (13) 

n 

with a suitable set of orthonormal basis functions 
{fn{k)}. The explicit form of the relevant functions used 
here can be found in the appendix. These can be trans¬ 
formed to the real space representation, 

= (14) 

^ k 

For a given function A^^(R, r), the coefficients an{Q) can 
be calculated from 

a„(Q) = ^ r)e-*‘3«/„(r). (15) 

R,r 

d. Parameters of the model - The model in Eq. 0 
contains a number of parameters, which we summarize 
here for clarity. The bare electronic structure is de¬ 
termined by the hoppings ti, t 2 and ts, and the filling 
n = 1/Ns cr ^i,cr by tbe chemical potential ji. Unless 
otherwise mentioned ti = 1 sets the energy scale. To 
get a rough idea about absolute scales we can think of a 
typical identification ti 300meV^ 3481K; however, we 
make no attempt for a quantitative theory in comparison 
with experiment here. The spin fluctuation spectrum has 


Sij = 2>g^x{ri 


(16) 


This purely static renormalization does not lead to any 
quasiparticle damping. However, it leads to a sizeable 
renormalization of the chemical potential and hopping 
parameters. Using the relative coordinate r = Vi — rj, we 
can relate the renormalized parameters {/i, U} to the bare 
ones by /i = /To + E(r = 0), ti = t? + E(|r| = 1) 

^2 = ^2 + ^ 3=^3 + E(|r| = 2), where 

t^ are the bare parameters. 

In our procedure we search for spontaneously ordered 
solutions of the self-consistency equations ^ and 0by 
initializing the calculations by a random field hij [see 
Eq. 01 - This is initially included in the Green’s func¬ 
tion in Eq. 0 and then set to zero from the second iter¬ 
ation on. We use a real space lattice with Ni = 32 sites 
in one direction and iteratively calculate the full Green’s 
function Gij and self-energy Some mixing of itera¬ 
tions and initial onsite homogenization is used to improve 
convergence. One example for a charge order solution is 
obtained with the parameter set t^ = 0.764, t^ = —0.33, 
4 = 0.154, and po = 0.393. Eor T = 0.05, F = ^“^ = 0.5, 
and g‘^ = 3, we obtain U « 1.0, t 2 ~ —0.32, ts ^ 0.11, 
p ~ —1.1, and the renormalized Fermi surface looks very 
similar to Fig. From the converged result for the 
expectation values Aij and Aare obtained as described 
in Sec. HH In Fig. § we plot the result for the nearest 
neighbor bonds in x— and ^—direction as obtained from 
A-^. This shows how an ordering pattern spontaneously 
appears and is stabilized after 30 iterations. 

The properties of the ordering pattern can either be 
identified directly from the real space representation in 
Fig-i or well understood by the decomposition into the 
basis function, Eq. (15). The largest coefficients are 
shown in Eig. Eor simplicity, we only show the co¬ 
efficients for a series of momenta Q = {Qx^Qy) iii fbe 
triangle, Qx G [0,7r], Qy G [0,Qx], labeled by uq. 
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FIG. 2: (Color online) Results for the static model. We plot the nearest neighbor bond values of Afj on a real space lattice 
initialized by a random field after 5 iterations (left) and spontaneously ordered results after 30 iterations (right). 



FIG. 3: (Color online) Plot of the largest an{Q) for different 
basis functions (n = 1,..., 13). The values of Q are labeled by 
ng, where nq = 25 corresponds to Q = Qo(l, 1), with Qo = 
3/8 in units ir/a. This value agrees well with the diagonal 
distance between hot spots. 


There is a clear maximum for as for nq = 25 which 
corresponds to Q = Qo{l,l), with Qo = 3/8 = 0.375 
in units ir/a. This value agrees well with the diagonal 
distance between the hot spots on the Fermi surface, and 
the ordering form factor is of the d-wave form, n = 3, 
fsik) = cos kx — cos ky (see Tablein the appendix). As 
mentioned we only show the result for a restricted set of 
Q vectors. The order in Fig. [^is really a superposition 
of ±(Qa,,Qy), ±(Qcc,—Qy) wave vectors. We conclude 
that our real space Eliashberg calculations are consistent 
with earlier work based on an instability analysisThe 
/c-space resolution for the finite size lattice is sufficient 
to resolve these features and the relation to the Fermi 
surface. The dominant instability is realized here as an 
ordered solution. 
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FIG. 4: (Color online) Plot of interacting Fermi surface in 
the dynamic model with weak coupling, ^ = 1.0. We plot the 
spectral density {pk{(^ = 0)) (left) and renormalized band 
structure (logpfc(cj) for clarity) for k along the trajectory 
F ^ M ^ A ^ F (right). The dashed line gives the non¬ 
interacting dispersion. The parameter values are given in the 
text. 



IV. RESULTS FOR DYNAMIC MODEL 


A. Homogeneous situation 


As discussed in the last section, in the static limit 
self-energy corrections only lead to a renormalization of 
the hoppings, but do not lead to any damping of quasi¬ 
particles. For a dynamic interaction such as in Eq. @ 
the effect is differentIn fact, this one of the ma¬ 
jor aspects of our work to include those effects on the 
single particle level. Scattering is particularly strong at 
the hot spots such that quasi-particle excitations can be¬ 
come suppressed. Here, we would like to investigate the 
effect of this on the charge order instabilities. The ba¬ 
sic questions are: How much is the diagonal order sup¬ 
pressed by the self-energy effects? Can the axial order 
become favorable over the diagonal one? In order to 
get some insights about the effects on the Eermi surface, 
we first consider the homogeneous situation and solve 
the Eliashberg equations self-consistently on the real fre¬ 
quency axis, Eq. (11). Eor clarity we first show results 


for weak interaction g = 1. The other model parameters 
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are t? = 1, = —0.32, = 0.11, /io = —1-1, ^sf = 0.5, 

ay^ = 1, r = 0.5, and T = 0.05. The result for the in¬ 
teracting Fermi surface, pk{oo = 0), for one part of the 
Brillouin zone is shown in Fig. (left). 

A broadening value r] = 0.08 was used to enhance clar¬ 
ity for the given /c-space resolution. Not surprisingly, 
the result resembles very much the non-interacting Fermi 
surface. On the right in Fig. we also display Pfc(cj) 
for k = {kx^ky) on the trajectory (0,0) ^ (7r,7r) ^ 
(0,7r) ^ (0,0) {r ^ M ^ X ^ F), which gives insights 
on the damping and renormalized band structure. We 
have added the non-interacting dispersion as a dashed 
line. Here again the result is almost identical with the 
non-interacting case. 

We also show the results for a stronger interacting case 
^ = 5 in Fig. where /i = — 1 and we show the temper¬ 
atures T = 0.1,0.05,0.02. Here, r] = 0.05 was chosen in 
the analytical continuation for the broadening. 

As we can see in Fig. the thermal part of the spin 
fluctuation spectrum can give considerable effects on 
the electronic spectrum and give features resemblant of 
Fermi arcs, which have been observed in experiments.^^ 
The effect is more pronounced for larger temperature 
as can also be deduced from Eq. 0- Hence, spec¬ 
tral weight is suppressed close to the hot spots where 
the scattering with the spin-fluctuations is particularly 
strong. Stronger arc-like features can be realized for 
larger values of the correlation length i.e., smaller val¬ 
ues of F. It is worth noting that non-selfconsistent cal¬ 
culations also give stronger arc features. However, self- 
consistent Eliashberg equati ons gen erically do not pro¬ 
duce real pseudogap features! ^^ * ^^ * ^^ * 

In Eig. we show for T = 0.05 the renormalized band 
structure. We find effects of Eermi velocity renormaliza¬ 
tion and redistribution of spectral weight. As we will 
discuss below this has consequences for the instability 
analysis and realization of ordered phases. 


B. Spontaneous order 

Within a similar procedure as described in Sec. Ing 
we have done extensive calculations to check for charge 
order solution for the model with the dynamic interac- 



FIG. 6: (Color online) Plot of renormalized band structure 
(log pfc(ct;) for clarity) at strong coupling g — 5.0. Other model 
parameter as described in the text. 


tion. The calculations were done for a real space lattice 
with Ni = 32 and a grid of Matsubara frequency large 
enough to capture relevant features. A random field hij 
is used for initialization. In the appendix, we describe 
some results which were obtained from calculations with 
field induced order. Eor most calculations we used the 
hopping parameters t? = 1, = —0.32, and tg = 0.11. 

By varying /io? we analyzed parameters with filling fac¬ 
tors n ^ 0.8 — 0.95 and temperatures down to T = 0.02 in 
units of The parameters of the spin fluctuation spec¬ 
trum where varied in a regime c^sf ^ 0.2 — 1, and F = 1/2 
and F = 1/4 was analyzed; we usually kept = 1. 

We also scanned a range of coupling strengths g. 

Generally, we found the appearance of the charge or¬ 
der to be strongly suppressed in a scheme with dynamic 
spin fluctuations as compared to the static case. One 
reason is the finite extent of the interaction in frequency 
space in contrast to the static case. Another reason is the 
single-particle renormalization effect of the self-energy. 
Charge order enters via an inhomogeneous modulation 
of in the self-consistent calculation, which ap¬ 
pears on top of generic variation of as function 

of Vi — Vj. This ordering tendency should be enhanced 
with the coupling g. However, the self-energy also has 
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FIG. 7: (Color online) Nearest neighbor bond values of Afj on real space lattice initialized by a random field after 5 iterations 
(left) and spontaneously ordered results after 45 iterations (right). Model parameter as described in the text. 


an effect to damp and renormalize single particle excita¬ 
tion, which leads to the opposite effect of a suppression of 
the order. This is a generic feature, which also appears 
for other Fermi surface instabilities, such as supercon¬ 
ductivity. As a consequence for many of the probed pa¬ 
rameters no charge ordered solution could be stabilized, 
even though some ordering features appear in intermedi¬ 
ate values of the iterations. For instance, for T = 0.05 
we did not find spontaneously ordered solutions for any 
of the parameters tested. 

In Fig. we present the result of a calculation at 
T = 0.02, where spontaneous charge order appears in 
the self-consistent Eliashberg equations. The param¬ 
eters are cjgf = 0-5, T = 0.5, with a coupling g = 
5. Since x{'^i ~ is fairly small this coupling 

strength is not as large as it might appear; for instance, 
the nearest neighbor self-energy reaches maximal values 
\ReT^{ri,ri -h x,iu;n)\ r\j 0.2 for these parameters. We 
show a result for /io = “I, which correspond to a filling 
of n 0.92 for the interacting system. The calculations 
are again initialized by a random field and the sponta¬ 
neous ordering pattern visible after 45 iterations is seen 
on the right in Fig.[^ The corresponding Fourier analysis 
is displayed in Fig.^ 

As in the static case we find a dominant d-form fac¬ 
tor, n = 3 and a diagonal wave vector with magnitude 
Q = Qq{1^1), with Qq = 3/8. There is also a minor 
component for n = 7, where f 7 {k) = cos{2kx) —cos{2ky). 
The wave vector matches the connection of the hot spots 
of the interacting Fermi surface well (as seen in Fig. [^. 
The diagonal order is dominant even though the spectral 
density is substantially suppressed near the hot spots. 
For different filling factors, we also find diagonal order, 
whenever a hot spot is present in the interacting Fermi 
surface, and the length of the wave vector changes ac¬ 
cordingly. 

In situations where there are no hot spots visible in the 
interacting Fermi surface, for instance, when the filling 



FIG. 8: (Golor online) Plot of the largest an{Q) for different 
basis functions (n = 1, • • •, 13). The values of Q are label by 
riQ, where uq = 25 corresponds to Q = (3o(l, 1)? with Qo = 
3/8 in units ir/a. This value agrees well with the diagonal 
distance between hot spots. 

is lowered such the Fermi surface becomes closed around 
(0,0), we found solutions with ordering tendencies with 
wave vector along the axis. The dominant form factor is 
n = 3, the d-form factor. For filling T = 0.02,n^0.8 the 
wave vector Q = ( 50 ( 150)5 with Qq = 1/2 can be found, 
which is very similar to the experimental wave length 
A ^ 4. The order is, however, quite small and not well 
established. It is possible that lower temperatures and 
additional ingredients are needed to stabilize the order 
well. 


V. DISCUSSION AND CONCLUSIONS 

Based on a model of fermions on a square lattice in¬ 
teracting with dynamic antiferromagnetic interactions. 
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we have analyzed general charge order solutions. We 
used self-consistent Eliashberg equations to compute self¬ 
energies and spectral functions. We showed that the 
spectral functions can be substantially modified in the 
interacting theory. In particular, structures resembling 
Fermi arcs can appear due to the strong scattering 
around the hot spots. However, no real pseudogap fea¬ 
tures are generated by self-consistent Eliashberg equa¬ 
tions. 

The appearance of charge order is studied based on 
a real space version of the Eliashberg equations, which 
are initialized by a random field. We only found spon¬ 
taneously ordered solutions at low temperature, e.^., 
T = 0.02ti, and for relatively large interaction strengths. 
Generically, the charge order posses a d-form factor and 
a diagonal wave vector related to the hot spots of the 
interacting Fermi surface. This result is in line with a 
number of previous stndies,^* ^^ * ^^ * ^^ * ^ ^^^ many of which 
however did not include dynamic interactions and self¬ 
energy effects, and studied instabilities from the normal 
state rather than ordered solutions. One of our main 
results is that the suppression of spectral weight at the 
hot spots by the dynamic model is not large enough to 
alter the behavior of the charge order. For large dopings, 
when the interacting Fermi surface has no hot spots, we 
found tendencies towards order with wave vector along 
the axis. This finding is in line with studies suppress¬ 
ing the hot spots, for instance, by pocket forma tion or 
including a phenomenological pseudogap SSHMMU 

In relation to the observation of charge order in the 
cuprate superconductors, we find a number of features 
consistent with the experiments. These include the dom¬ 
inant d-form factor and the behavior of Q with dop¬ 
ing. The direction of the wave vector is, however, not 
correctly reproduced. This is likely related to the fact 
that the self-consistent Eliashberg theory does not in¬ 
clude strong enough pseudogap features, and therefore 
still supports an instability related to the hot spot. As 
such, the present approach does not provide a compre¬ 
hensive explanation for the experimental findings. Never¬ 
theless, with suitable modifications our real space Eliash¬ 
berg approach might useful for future applications, as it 
can easily be extended to include superconductivity, mag¬ 
netic field effects, impurity effects, coupling to phonons 
and spin order. 
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n 

fn{k) 

1 

1 

2 

COS kx + cos ky 

3 

cos kx — cos ky 

4 

sin kx + sin ky 

5 

sin kx — sin ky 

6 

COs(2/Ca;) -h COs{2ky) 

7 

cos(2fc) — cos{2ky) 

8 

sin(2/cx) + sin(2A:2/) 

9 

sin(2/ca;) — sm(2ky) 

10 

2 cos kx cos ky 

11 

2 sin kx sin ky 

12 

2 sin kx cos ky 

13 

2 cos kx sin ky 


TABLE I: Relevant basis functions 


Appendix 
1. Basis functions 

Here we collect the basis functions which were used for 
the charge order analysis. We focused on the 13 basis 
functions fn{k) as shown in TableThe corresponding 
function of the real space representation fn{r) can be 
easily calculated analytically from Eq. ( p^ . They include 
form factors up to next-nearest neighbors. Higher order 
functions could be in principle be considered but they do 
not play a role in our analysis. These functions fn{r) are 
used to compute the coefficients an{Q) in Eq. ( p!^ . 


2. Field induced order 

Here, we briefly summarize some results from a com¬ 
plementary analysis to Sec. IV. We studied charge or¬ 
der when a small finite field is present in the self- 
consistent calculations. We focused on the situation with 
d-form factor and assumed 03 (Q) = 03 (—Q). We chose 
03 (Q) = 0.001 keeping the field finite in 0 during the 
self-consistent calculation and scanned over various vec¬ 
tors Q. A response function characteristic of the suscep- 
tility to charge order can be defined as x(n, Q) = j . 

First of all we can consider the situation for = 0, such 
that T^ij = 0. Glearly there cannot be spontaneous order 
in this situation, however, the susceptibility for different 
wave vectors can be analyzed. We find that the field 
induced order as measured by X^(3,Q) is strongest for 
Qd = {Qo^Qo)^ where Qo approximately connects the 
hot spots, and there is also a local maximum along the 
Qx axis, for Qa = (Qo, 0)- Within the real space calcula¬ 
tions with Ni = 32 one finds x^(3, Qa)/x^(3, Qd) — 0-78. 
Hence the diagonal order is favored. These results are 
clearly expected from the instability analysis based on 
unrenormalized fermion^^ where these quantities can be 









computed from the corresponding fermionic bubbles and 
form factors. 

For finite g the self-energy effects play an important 
role. As discussed the ordering tendency is induced 
form an inhomogeneous contribution of the real part of 
In contrast to the static case T>ij{iuJn) 0 for 
large iujn so this is only a contribution at small ujn • How¬ 
ever, the self-energy also produces other effects, such as 
renormalization of the band and suppression and shift 
of spectral weight as already discussed in the section 
on the homogeneous calculations. These effects mostly 
lead to a reduction of the ordering tendency. This is in 


fact common for weak coupling instabilities which are 
strongest without self-energy corrections. Hence, the 
combination of these effects means that the field induced 
ordering susceptility is not necessarily enhanced for finite 
g. In our analysis for different parameters and interac¬ 
tions we found that typically x(3, Qd) with the diago¬ 
nal wave-vector is largest and can be enhanced over the 
non-interacting value. For the parameters studied diag¬ 
onal response was found to be larger than the x(3, Qa)i 
i.e., instabilities along the axial direction, consistent with 
what has been discussed in the main text. 
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